Conductance of a molecule with a center of mass motion 
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We calculate the zero temperature conductance and characteristic correlation functions of a molecule with a 
center of mass (CM) motion which modulates couplings to the leads. In the first model studied, the CM vibra- 
tional mode is simultaneously coupled to the electron density on the molecule. The conductance is suppressed 
in regimes corresponding to non-integer occupancy of the molecule. In the second model, where the CM mode 
is not directly coupled to the electron density, the suppression of conductance is related to the dynamic breaking 
of the inversion symmetry. 
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I. INTRODUCTION 

Interest in systems with electron-phonon interaction has 
recently reemerged, as a consequence of advances in ex- 
perimental techniques enabling the investigation of the 
subtle interplay of electronic and vibrational degrees of 
freedom by means of measuring the conductance of the 
moleculesi2iiiiSi^ii2i2ii2iiiiiii^. Consequently, a considerable 
amount of work was done in pursuit of understanding the 
equilibriumiiii^ii^iii and non-equilibrium. '^ '^-^^ properties of 
these systems. The analysis of low-temperature properties 
of such systems is challenging as strong repulsive interaction 
among electrons confined to molecular orbitals leads to sur- 
prising results such as the interplay of the Kondo physics and 
molecular vibrations, which are intractable to conventional 
perturbational approaches. Investigation of the molecules 
which, when attached to the leads, may posses both internal 
and center-of-mass vibrational modes, makes the calculation 
of various properties even more involved. 

Here we consider a molecule undergoing a linear transport 
measurement - a system which consists of a molecule and at- 
tached leads. The shape of the molecule and its position with 
respect to the leads oscillate (the setup is sketched in Fig^. 
Hence the orbital level of the molecule - for transport through 
molecules it is sufficient to consider a single orbital due to 
wide inter-level spacings^ - and the tunneling amplitude are 
both modulated. Here we treat two possibilities: (i) in case 
I (proposed in Ref. I21I) the vibrational mode that modulates 
the tunneling amplitude is also coupled to the electron den- 
sity at the molecule; (ii) in case II (treated also in Ref. l22lh 
the molecule possesses in addition to the CM mode a breath- 
ing mode which modulates the energy level of the molecular 
orbital. 

Experimentally, CM modes are discerned in nonlinear 
transport measurements as side-bands in conductance which 
correspond to frequencies which do not match any of the 
eigenmodes of a given molecule. While the relevance of the 
separate couplings to two vibrational modes in model II is in- 
tuitively clear, the relevance of the simultaneous modulation 
of tunneling and coupling to the electron density in model I 
is less clear and deserves further exploration. We note that 
the model I could be relevant for the nonlinear transport mea- 
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Figure 1: (Color online) Schematic plot of the model devices. In 
model I the CM vibrational mode which modulates the tunneling am- 
plitude t' is coupled also to the local charge density n. In model II the 
CM displacement modulates only the tunneling amplitude whereas 
another (breathing) vibrational mode is coupled to the local charge 
density. 



surement where the electric field would directly couple to the 
electron density at the molecule, or in the presence of some 
impurities with residual electrical fields even in the (near) lin- 
ear transport measurements. 

The case of a molecule without the CM oscillations de- 
scribed by the Anderson-Holstein Hamiltonian has already 
been analyzed by several authors (see, e.g. Ref. [i3 and ref- 
erences therein). The dominant effect of the coupling to the 
breathing mode is the reduction of the repulsive interaction 
among electron pairs occupying the molecular orbital. When 
the reduced repulsion J/off is positive, a gate-voltage sweep re- 
veals the area (of width of order Ucs) of enhanced linear con- 
ductance due to the Kondo effect. With increased electron- 
phonon coupling and thereby decreased Ucs this area is di- 
minished. In this paper we show that the CM oscillations sup- 
press the conductance in a new way which cannot be described 
by an effective model with reduced repulsion. 

This paper is organised as follows. Section II introduces 
two related electron-phonon models. In Sec. Ill numerical 
methods based on the variational wave function are presented 
and some symmetry properties of the models are considered. 
In Sec. IV results for linear conductance for both models are 
presented and some ground state properties of the molecular 
system are explained. Results are summarised in Sec. V. 



2 



II. MODELS 

Hamiltonians for both models consist of the electron, the 
phonon and the electron-phonon coupling terms, H = He + 
Hp + Hc-p. In the electronic part of the Hamiltonian, 

He = Hi, + Hr + edUd + Und]ndi - t'{vL + vr), (1) 

^L.R = —t J2ieL R "^I'-i+i + ^-^^ describe the noninteracting 
tight-binding left and right leads, respectively, and is the 
energy level of the molecular orbital occupied by n = (n^) 
electrons where nd — J2a ''^d.cr with Uda — d^d^. The char- 
acteristic scale of the repulsive interaction is U and the tun- 
neling part with characteristic tunneling rate t' consists of op- 
erators vi — cjj^do- + h.c. for / = L, R which couple 
the molecular orbital to the first site in the left and right leads, 
respectively. In He only the coupling between the molecular 
orbital and even combinations of the left and the right lead 
orbitals (i.e., the symmetric channel) is present. 

Vibrational modes are described with H^ = Vtia\ai and 
H^ — ilifljai + Vt2a\a2 for models I and II, respectively, 

and a\ ^ are the phonon creation operators. The corresponding 
electron-phonon interaction is given by 

I : H\_p = t'giv-L - vh)xi + Mind - l)xi 
II: Hllp=t'gidi^-dn)xi+M{nd-l)x2, (2) 

where xi^2 = a! 2 + '^1.2 are the displacement operators. 



III. METHODS 

The ground state properties are determined using the Gun- 
narsson and Schonhammer projection-operator methodii2i2i. 
Here the Hamiltonian is diagonalized in the basis 



l*A{™„}) = P> 



|5). 



(3) 



which consists of projectors P\ , Pq — (1 — 71^1) (1 — ndi). 
Pi = J2a^d.a (1 - nda), and P2 = nd-\ndi and additional 
operators involving the operators in leads (for example, P3 = 
PqVlPi), which are applied to the state |0) corresponding to 
the phonon vacuum and the ground state of the auxiliary non- 
interacting Hamiltonian H = idnd + ^iVl + t'^VR + Hi, + Hr 
with renormalized local energy id and hopping parameters 

j^. In the model I a = 1 and a — 1,2 for the model II. 
Due to the CM displacement which induces the coupling to 
the asymmetric channel, the renormalized couplings to the left 
and right lead are not necessarily equal. The chemical poten- 
tial was set to the middle of the band. A sufficient number of 
variational operators (up to ~ 40) and excited phonon states 
(up to ~ 40) were used in order to obtain converged results in 
the parameter regimes presented here. 

Zero temperature conductance was calculated by two re- 
lated methods, both based on the ground state wave func- 
tion expressed in the basis Eq. Q. In the first method 



it was calculated from the sine formula (SF)^^'^^, G = 
Gosm^[{E+ - E^)N/4t], where Go = 2e^/h and E± are 
the ground state energies of a large auxiliary ring consisted of 
N non-interacting sites and an embedded interacting system 
(molecule), with periodic and anti-periodic boundary condi- 
tions, respectively. In the second method we use the effec- 
tive parameters i'j^R,£d of the auxiliary effective noninteract- 
ing Hamiltonian H and the conductance is then given from 
the corresponding Green's function (GF)24. The advantage of 
the GF method is that the corresponding GF can be evaluated 
for infinite leads so the difficulties with convergence of SF 
with number of sites are avoided. However, the SF method 
is robust and it depends only on the accuracy of the ground- 
state energy which improves with the size of the basis in a 
transparent way. By comparing results of both methods we 
checked for the consistency and the convergence. From the 
effective parameters the conductance can also be calculated 
using the Friedel sum rule valid in the limit of large band- 
width {t :§> U, M, rii, ...). When the coupling of the system 
between left and right lead differ, conveniently parameterized 



as tj 



1 + b)t', = (1 — b)t', the conductance is given by 



G_ _ jl-b^y 



2 sm —n. 



(4) 



Note that conductance may be less than unity when both the 
symmetric and the antisymmetric channel participate in the 
transport {i.e., when & ^ A 6 7^ ±cxd). 

The asymmetry stems from the coupling to the CM dis- 
placement, hence the expectation values of the CM displace- 
ment x — {xi), and the calculated asymmetry factor b are re- 
lated. More precisely, the following exacU^iS relations hold: 

/ ^ 2,gt' _ _ . 2M , 

II : (xi) = - ^^r), {x2) = -^{n - 1). (5) 

ill ^'2 

When n ~ 1 the displacement x is proportional to the activity 
of the antisymmetric channel (ul — vr) for both models. 

Some general conclusions may already be drawn from con- 
sidering the particle-hole symmetry. Model I is invariant with 
respect to the particle-hole transformation cj^. cia ,d^^d 
with (e + U/2) ^ -(e + U/2), t -t, t' -> -i' and 
xi —xi. Hence at the point of particle-hole symmetry 
(where n = 1), x = 0, therefore 6 = and conductance is 
unity according to the Friedel sum rule Eq. (|4}. The particle- 
hole transformation of model II yields X2 ~X2 with no 
bounds on xi. The conductance of model II in the symmet- 
ric point therefore cannot be deduced from considering the 
particle-hole symmetry alone and may be less than unity there. 



IV. NUMERICAL RESULTS 

Throughout this paper we analyze the system with positive 
Ues ^ U ~ 2IVP/na and show the results for U = IGF, 
M = ria = 2.5r, C/off = 5r,r/t = 2t''^/t^ = O.OS, where 



^l(b)' I— „Jl 
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Figure 2: (Color online) (a) Conductance as calculated from the 
Green's function corresponding to renormalised parameters (GF, full 
line), sine formula method with = 3000 (squares), and Friedel 
sum rule (circles) for molecules described by model I with g = 0.4. 
For comparison, conductance for a molecule without coupling to the 
CM mode, g = (dashed), (b) Local charge occupancy n — 1 
(full), the displacement x/2 (short-dashed) the asymmetry param- 
eter b (dotted) for a molecule described by model I with g = 0.4, 
and the occupancy n — 1 for <? = (long-dashed). 

a — 1,2 for model I and II, respectively. The scales of the 
values used correspond to typical molecular transistorsi. 

A. Model I 

In Fig. 13 a) we show that all methods yield the same re- 
sult for conductance within the model I (the same is true for 
model II discussed later). The area of enhanced conductance 
is diminished in comparison to the g — Q case. In general, 
the conductance curves calculated using different methods are 
very close. However, due to the finite bandwidth used here 
minor discrepancies in results applying the Friedel sum rule, 
Eq. 0, might occur when compared to other methods in the 
charge transfer and empty orbital regimes, Ed + U/2 > t. 

The results for the asymmetry factor b, the occupancy n — 1 
and the displacement x are shown Fig. |2jb). As discussed 
above, the displacement at the symmetric point has to van- 
ish. Correspondingly, the asymmetry 6 = there. Near the 
particle-hole symmetric point we find the relation between the 
displacement and the asymmetry factor to be approximately 
linear. 

We now systematically explain the properties of the ground 
state for model I starting from a Hamiltonian with g — 0. In 
Fig-El^) we plot the conductance. Surprisingly, the influence 
of increasing g on the conductance is similar to simple reduc- 
tion of U as is the case with the Anderson-Holstein model 
where the electron-phonon coupling term M effectively re- 
duces U. However, as shown in Fig. |3jb) the occupancy n 
does not exhibit the sharp transition from n ^ 2 to n ~ 
as seen in the Anderson model with reduced Ucs, where elec- 




Figure 3: (Color online) (a) Conductance of a molecule described by 
model I for various CM couplings g and fixed Holstein coupling A/, 
2AfVfii = U/2, U = lOr, M = 0,1 = 2. sr. (b) The occupancy 
n— 1 for various (;. Forg^ = 0.25 also the average displacement a;/2 
(dashed-double-dotted) plotted, (c) Fluctuations of charge An^ = 
{(rid — n)^) and (d) displacement Ax^ = {{xi — x)"^). 



tron (or hole) pairs are preferred in the ground state (see e.g. 
Ref.Djv and the transport is dominated by the pair tunneling^it 
In the strong coupling regime the occupancy preferentially 
takes on values near half integer instead. This is because to 
the lowest order only the average displacement couples the 
electronic and phononic part, and the average displacement 
vanishes when n = 1. Also, for states with well defined occu- 
pancy n ~ 0, 1, 2, the hopping matrix element has a vanishing 
weight. 

To illustrate this further we determined the charge fluctu- 
ations on the molecule. In general the range coiTesponding 
to increased charge fluctuations An^ for 5 > [Fig.|3c)] is 
extended compared to the g = case. For large g the charge 
fluctuations still exhibit a minimum at the particle hole sym- 
metric point, in contrast to the negative-?/ (corresponding to 
strong coupling to the breathing mode) case where the charge 
analog of the Kondo effect with increased charge fluctuations 
evolves'^. The straightforward application of the fluctuation- 
dissipation theorem relates A-n? ^ —dn/ded, which does not 
apply to this system in some parameter regimes. On the other 
hand, the displacement fluctuations shown Fig. |3jd) are in- 
creased for large g at the symmetric point as in the Anderson- 
Holstein model 

In Fig.|3jb) we also plot the average displacement x/2 of 
the molecule for g = 0.5. According to the relation in Eq. (|5Jl, 
the difference between the occupancy (times 2Af/17i) and the 
displacement x is determined by the difference between ex- 
pectation values of the tunneling into left and right leads. At 
the symmetric point, this difference vanishes. 

To summarize, in the particle-hole symmetric point the con- 
ductance is unity, contrary to the corresponding results pre- 
sented in Ref. |2lj where a dip in conductance was reported 
and related to the non-applicability of the Friedel sum rule 
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Figure 4: (Color online) (a) Conductance of a molecule described by 
model II for g < Qc (dashed) and Qc < g < g'c (full). Sli = 0.2t. 
Inset: The occupancy for the two cases. For g > gc also the CM 
displacement x is shown (dotted). For g < gc the CM displacement 
vanishes, a:: = 0. (b) Fluctuations of charge and of CM displacement. 
The curve of charge fluctuations for g = 0.5 is valid for all g < g^. 
(c) The conductance in the symmetric point as a function of g. 



due to the breakdown of the Fermi liquid. Our results con- 
firm that the conductance can be expressed with the Friedel 
sum rule, if left-right asymmetry is correctly taken into ac- 
count, Eq. (0}. The displacement of the oscillator in model I 
is determined by the occupancy. The displacement increases 
coupling to the anti-symmetric channel and thereby the charge 
fluctuations which destroy the Kondo correlations with well 
defined charge and enhanced conductance. Further, the dis- 
placement gives rise to asymmetry, which means conductance 
is additionally suppressed for general ^ —U/2. The width 
of the conductance peak is seen to decrease approximately as 
U - 2AP /VLi - /(Af, fii, t')g^ until this quantity is reduced 
to zero. 



B. Model II 

We now turn to the results within model II. For small 
enough values (to be quantified below) of the CM coupling 
g < Qc the conductance is not significantly affected as unlike 
in model I the asymmetry is not determined (directly) by the 
occupancy of the molecule. The average displacement van- 
ishes, correspondingly there is no asymmetry in the effective 
model; the width as well as the shape of the conductance peak 
and also the correlations exhibit only minor renormalisation 
compared to the g = case. Conversely, for g > gc a soft 
mode emerges and the model has a transition reminiscent of 
dynamical Jahn-Teller distortion'^ as within a simplified treat- 
ment of the related model^^. This is signaled in our approach 
with non-vanishing x and asymmetry b in the effective model 
H and correspondingly with qualitatively modified correla- 
tions. For this regime additional variational basis functions 
might be required^". For CM coupling large enough, g > g'^. 



only the asymmetric channel is active. This is a consequence 
of the renormalization: only the channel with the strongest 
Kondo coupling remains active at low temperatures, in accord 
with Ref. 22, and with coiTesponding unity conductance. 

In Fig. 13(a) we plot the conductance and in Fig.|3(b) the 
fluctuations of occupancy and the displacement as a function 
of ed for two distinct gc < g < g'c and g < gc and fixed 
ili/t — 0.2. Whereas for g < gc the conductance and cor- 
relations are almost unaffected by the CM coupling (compare 
with the g — curves in Figs. |2 0), the g > gc result is 
markedly different. In the latter, the conductance is less than 
unity at the symmetric point and is zero at the point coiTe- 
sponding to 5 = 1. In the inset of Fig. |4] (b) we also plot 
the occupancy and for g > gc the CM displacement. For 
large enough, the CM displacement (and correspondingly the 
asymmetry) is zero. The precise value for gc hence varies with 
the gate voltage and is related to the charge fluctuations of the 
bare model: increased charge fluctuations induce a transition 
to the state of spontaneously broken symmetry. 

In Fig. 13(c) we plot the conductance as a function of cou- 
pling to the CM mode. For small g the conductance is unity 
as only the even channel is active, 6 = 0. By increasing g, 
the conductance drops sharply at gc which is found to be set 
by the relation gIt/Vli = c(t'), where c(t') '-^ 2.5 for the 
results presented here^. Note labels A and B indicating the 
values of g for which the results are shown in Figs. |4la,b). 
For large fii both channels are active with coiTesponding 
suppressed conductance only for a narrow range of g > gc 
while for progressively decreasing fli the range of g with sup- 
pressed conductance widens. As in typical devices considered 
in experiments^ fli < T < U the anomalies in conductance 
of the type considered here would would not be sensitive to 
the precise value of the coupling to the CM mode provided 
that it is large enough. 



V. CONCLUSIONS 

In conclusion, we analyzed the properties of molecules with 
a center-of-mass vibrational mode in linear transport mea- 
surements described by two distinct Hamiltonians. In the 
first case, where the coupling to the CM mode simultane- 
ously modulates the molecular energy level, the result of the 
coupling to the CM mode is a reduced range of gate volt- 
ages with enhanced conductance. As a direct consequence of 
the particle-hole symmetry the conductance in the symmetric 
point remains unity. 

In the second case where the coupling to the CM mode 
modulates only the tunneling matrix elements and additional 
breathing mode couples to the electron density, the CM mode 
is of minor relevance for the transport until a critical value of 
this coupling gc is reached. For values of the coupling near 
the critical value anomalies due to the development of the soft 
mode arise signaled by the conductance which is less than 
unity. The conductance for these values is within the effec- 
tive model suppressed for a wide range of gate voltages. 

We explained the relation of CM displacement to the trans- 
port of electrons for both models and substantiated the expla- 
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nation by the numerical data for occupancy, occupancy fluc- 
tuations and also related phonon expectation values. 

Finally, we may make some remarks regarding the experi- 
mental relevance of results presented here. For linear regime 
without residual electric fields it does not seem likely that 
the model I would provide an accurate description of the ex- 
perimental situation. By analyzing model II which applies 
then, we have shown that for small couplings to the CM mode 
there are no discernible effects on the linear conductance (we 
should note here that this statement would not be changed if 
additional breathing modes were introduced). At a sufficient 
coupling strength this picture changes: the displacement of 
the system from the point of inversion symmetry becomes in- 
creasingly sensitive to external perturbations which would re- 
sult in a displacement of molecule and strong suppression of 
conductance. 

To check for this behaviour in experiments would require 
the following steps. First devices in which CM side-peaks 
are seen would have to be produced, possibly by using dif- 



ferent molecules, which would give in other respects identical 
devices besides having different frequencies of the CM mo- 
tion. The fingerprint of the physics considered here would 
then be the suppression of conductance which would predom- 
inantly occur for the devices with lower frequencies of the CM 
motion. As this procedure would require near-atomic similar 
lead configuration near the molecule the systematic experi- 
mental confirmation of the analyzed effects seems currently 
hard to achieve. We speculate nevertheless that this effects 
might have contributed to the anomalously suppressed con- 
ductance in some of the devices considered in, e.g., Ref.^S, 
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